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ABSTRACT 

Aims. We investigate the star formation ocurring in the region towards IRAS 07141-0920 contained in the HII region Sh2-294 (S294). 
We report the discovery** of a new young double stellar cluster, and we describe its properties. 

Methods. High-resolution optical UBVRI and Ha images obtained with ALFOSC mounted on the Nordic Optical Telescope (NOT), 
near-infrared JHKs images obtained with NOTCam at the Nordic Optical Telescope, and VLT/ISAAC images obtained through the 
H 2 (2.12 yum) filter were used to make photometric and morphological studies of the point sources and the nebula seen towards Sh2- 
294. 

Results. The optical images reveal an emission nebula with very rich morphological details. The nebula is composed mainly of ultra- 
violet scattered light and of Ha emission with local morphological features conveying an impression of great dynamical complexity. 
Contrasting with the bright parts of the nebula, opaque patches are seen. They appear elongated, filamentary, or clumpy. Our optical 
photometry confirms that the identification of the illuminator of the nebula is likely to be a B0.5V star located at a distance of about 
3.2 kpc. Our high-resolution near-IR images reveal an embedded cluster, extending for about 2 pc and exhibiting sub-clustering: a 
denser, more condensed, sub-cluster surrounding the optical high-mass B0.5V illuminator star; and a more embedded, optically in- 
visible, sub-cluster located towards the eastern, dark part of the nebula and including the luminous MSX source G224. 1 880+0 1 .2407, 
a massive protostellar candidate that could be the origin of jets and extended features seen at 2.12 jim. The double cluster appears to 
be clearing the remaining molecular material of the parent cloud, creating patches of lower extinction and allowing some of the least 
reddened members to be detected in the optical images. We find 12 main-sequence (MS) and 143 pre-main sequence(PMS) members 
using three different methods: comparison with isochrones in optical colour-magnitude diagrams, detection of near-IR excess, and 
presence of Ho- emission. The most massive star fits a 4 Myr post-MS isochrone. The age of the optically selected PMS population is 
estimated to be 7-8 Myr. The IR-excess population shows sub-clustering on scales as small as 0.23 pc and is probably much younger. 

Key words. Stars: Formation - ISM: Clouds - ISM: Individual (IRAS 07141-0920) - ISM: Individual: Sh2-294 (S294) - Infrared: 
stars - (ISM:) dust, extinction - Stars: pre-main sequence 



1. Introduction 

Galactic nebulae are often times clear signposts of regions of 
star formation due to the proximity of the young stars to the par- 
ent material of the molecular cloud that originated the stars. The 
presence of nebulosity has been used to find and identify star for- 
mation regions (e.g. |Yun & Clemens |1994| |Tapia et al. | 1997 ; 



Rubio et al. 1 1998 1. In addition, young stars are the source of en- 



* Based on observations made with the Nordic Optical Telescope, op- 
erated on the island of La Palma jointly by Denmark, Finland, Iceland, 
Norway, and Sweden, in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. Based on obser- 
vations collected at the ESO 8.2-m VLT-UT1 Antu telescope (program 
66.C-0015A). 

** A new publication (Samal et al. 2007, see references in the text) 
with a similar analysis of a field in this region appeared when our paper 
was in the refereeing process. 



ergetic winds, jets or outflows, whose interaction with their sur- 
rounding cloud material produces shocks, which are denounced 
by nebulae with characteristic shapes (e.g. bow shocks, HH ob- 
jects (e.g. Reipurth & Bally |2001|l), and ultimately leads to the 



clearing of the cloud material and emergence of the young stars 
out of the cloud core. 

The morphology of nebulae seen towards star forming re- 
gions can reveal the processes that have been ocurring in the 
interior of the parent cloud. Diffuse reflection nebulae, both in 
the optical and in the near-infrared are the result of illumination 
from bright sources more or less embedded in the cloud. Optical 
emission nebulae corresponding to excited gas can be detected 
using narrow-line filters, the most common one being Ha (e.g. 
|Fich, Dahl, & Treffers |1990) i. 
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Small groups or isolated young low-mass stars have limited 
dynamical effect on their surroundings. On the other hand, rich 
clusters of young stars or high-mass young stars will often pro- 
duce spectacular nebulae. Massive star formation is currently 
being addressed using a variety of the new instruments made 
available in the last years. Unlike the low-mass family of young 
stars, the relatively short timescale for their accretion and pre- 
main sequence evolution implies that massive stars reach the 
main-sequence still embedded in their parent molecular cloud 
cores. On the other hand, their strong hard radiation ionizes 
their surroundings creating a ultra-compact H II (UCHII) re- 
gion, which may represent the earliest recognizable stage in the 
life of massive stars ( jChurchweli 2002). UCHII regions have 



been identified using Ha surveys (e.g. Rodgers, Campbell, & 
Whiteoak 1960 1, radio surveys, and looking for characteristic 



IRAS colours (e.g. [Wood & Churchwell |19 89). More recently, 
they have been studied using a combination of near and ther- 
mal infrared surveys and K-band spectroscopy ( |Bik at aL~| 2003, 
|2006) . 

IRAS 07141-0920 is an IRAS source with colours close to 



those of ultra-compact HII regions ( Wood & Churchwell 1989 ). 



It is seen towards a large, extended optical nebula. Rodgers. 
Campbell, & Whiteoak |(T%0|) detected Ha emission, and Fich. 



Dahl, & Treffers | ( |1990> detected the Ha line at v(LSR) of 33.5 
km s" 1 . The source has been first cataloged as an HII region by 



Sharpless ( 1959] > in his catalogue of HII regions (source Sh2 



294 or S294). Using UBV p hotoelectric photometry, Moffat, 



Jackson, & Fitzgerald ( 1979} derived a photometric spectral 



type of B0.5V for the hottest star in this HII region, which they 
assume to be the ionizing source. Their estimated distance was 
of 4.6 kpc. Extended radio continuum emission from S294 at 
1 .46 GHz has been detected and mapped, resulting in a diameter 
for the HII region of at least 7x7 arcmin 2 ( |Fich |l993| . A strong 
CO 7=1-0 molecular line transition towards IRAS 07141- 
0920 h as also been detected, at a radial velocity v(LSR) of 32 
km s" 1 dBlitz, Fich, & Stark 



1982 



Wouterloot & Brand 1989) 



They derived a kinematic distance of 3 kpc, much smaller than 
the value of 4.6 + 1 .5 kpc derived by optical photoelectric pho- 
tometry (|Moffat, Jackson, & Fitzgerald ||1979) . |Wouterloot &| 



Brand (1989 ) also note that the CO line is flat-topped (non- 



gaussian), a hint that the molecular gas is not in a quiescent 
state, possibly with turbulence or other line broadening mech- 
anism present. Because they have a single observation and have 
not mapped the cloud, there is no estimate of the extent of 
this molecular cloud seen towards IRAS 07141-0920. Searches 
for masers resulted in non-detections (water masers: Henkel. 



Guesten, & Haschick |(fl 986); methanol mas ers: MacLeod et al. 
(1998 ), Codel la et al. | ( |1995) ). |Bica et al. |p003) reported the 
presence of a 2MASS infrared, loose star cluster towards Sh2- 
294. |Samal et al. | (| 2007| l provides a multivalength analysis of 



the Sh2-294 region describing the distribution and nature of the 
young stellar objects, and the morphology of the thermal ionized 
gas and dust emissions^] Here we present an optical and near- 
infrared study of the extended optical nebula and the embedded 
stellar population seen towards Sh2-294 (S294) / IRAS 07141- 
0920. We confirm the presence of a young embedded cluster 
which appears to be actively in the process of destroying its crad- 
dle, that is dissipating the surrounding parent molecular cloud. 
Most cluster members are still embedded and are either not seen 
or are very faint on our deep optical CCD images. 



1 This article was published while our manuscript was in the referee- 
ing process. 



Section 2 describes the observations and data reduction. 
In Sect. 3, we describe the properties of the extended nebula. 
Section 4 describes the stellar density distribution and the pres- 
ence of a young embedded cluster. In Sect. 5, we present our 
photometric results, estimate the reddening, distance, and age 
of the cluster, and make cluster membership assignments. In 
Section 6, we discuss the spatial distribution of the cluster mem- 
bers and present a possible star formation scenario for this re- 
gion. A summary is given in section 7. 



2. Observations and data reduction 

2.1. Optical observations 

Optical (UBVRI and Ha) observations were carried out dur- 
ing the night of 2006 November 29 with the NOT (Nordic 
Optical Telescope), operating at the ORM (La Palma) under very 
good seeing conditions (0.6 arcsec). The ALFOSC (Andalucia 
Faint Object Spectrograph and Camera) CCD camera in imag- 
ing mode was used (with a plate scale of 0. 19 arcsec/pixel). Short 
and long exposure images were obtained through each band: 
short exposures with integration times of 10, 10, 5, 5, 5, and 
60s for UBVRI and Ha (FWHM=3.3 nm) bands, respectively; 
and long exposures with integration times of 600, 300, 300, 300, 
300, and 1200 s, for the same bands, respectively. 

The observations were reduced with the adequate routines 
in the IRAF^] package. In each of the UBVRI -bands, three long 
exposures were co-added into a single image. Aperture and PSF 
photometry was performed in both the short and the long expo- 
sure images. The resulting magnitudes were then joined together 
to produce a final list of instrumental magnitudes. 

Frames of a standard field were also secured at different air- 
masses. The field includes 12 stars in the SA98, around star 652 
( |Landolt |1992| l. These observations were used to transform the 
instrumental magnitudes to the standard system, with the pro- 
cedure described by Delgado et al. ( 1998), based upon the as- 



sumption of multilinear dependence of standard magnitudes on 
airmass, instrumental magnitude, and an appropriate colour in- 
dex. 



2.2. Near-infrared observations 
2.2.1. NOTCam JHK S observations 

Near-infrared (J, H and Ks) images were obtained on 
2007 January 29 using the Nordic Optical Telescope near-IR 
Camera/spectrograph (NOTCam^j The detector was the 1024 
x 1024 x 18 micron Hawaii engineering grade array. The wide 
field camera (0.234"/pix) was used, and the observations were 
performed with a ramp-sampling readout mode. Every sky posi- 
tion was integrated for 30 sec, reading the array every 5 sec and 
using the linear regression result of these 6 readouts. For each 
filter, 9 dithered sky positions were observed. Because of vary- 
ing sky conditions, only the best frames were selected, resulting 
in total on-source integration times of 210, 180, and 150 seconds 
for J, H, and K s , respectively. 



2 The Image Reduction and Analysis Facility (IRAF) is distributed 
by the national Optical Astronomical Observatory, which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Science 
Foundation. 

3 See URL http://www.not.iac.es:/instruments/NOTCam/ for more 
information on NOTCam. 
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The images were reduced using a set of own NOTCam 
scripts to correct for bad pixels, subtract the sky background, 
and flatfield the images. Differential twilight flats taken a few 
nights apart were used to correct for the pixel to pixel variations 
of the response. The selected images were then aligned, shifted, 
trimmed, and median combined to produce a final JHKs image 
of 3.7' x 3.4' size. 

Point sources were extracted using daof ind with a detection 
threshold of 6<x, and a few additional sources were added in by 
hand. Aperture photometry was made with a small aperture (ra- 
dius = 4 pix, which is about the measured FWHM of the point 
spread function) and aperture corrections, found from 5 bright 
and isolated stars in each image, were used to correct for the 
flux lost in the wings of the PSF. The error in the determination 
of the aperture correction was < 0.02 mag in all cases. These er- 
rors are added to the value MERR which is output by the IRAF 
task phot. A total of 491, 496, and 494 sources were found to 
have fluxes in J, H, and K$ , respectively, and errors cr Ks < 0.35 
mag in the Ks band (the noisiest and least sensitive). All these 
are included in the table of photometric results for complete- 
ness (Table 1, see below), but we include only the 444 NOTCam 
sources with o~ks < 0.2 mag in our analysis in this paper. 

The brightest star in the field S294-B0.5 (ID 717) had peak 
pixel values outside the linear range of the detector both in the 
H and K$ bands. We therefore used only the first of the 6 ramp- 
sampling readouts, with an effective integration time of 5 sec- 
onds for this source. 

We hav e used the 2MASS Al l -Sky Release Poin t Source 
Catalogue ( |Skrutskie et al. ||2006] |Cutri et al.~||2003[ ) to cali- 
brate our observations both in terms of photometry and astrom- 
etry. The plate solution was found using 27 isolated and bright 
2MASS stars over the FOV and the iraf task ccmap with the 
"tnx" projection allowing for geometric distortion and a 3rd or- 
der polynomial with full cross terms. NOTCam XY positions 
were then transformed to RA and DEC, and all 2MASS sources 
were matched to all NOTCam sources using a search radius of 
0.36", which corresponds to three times the positional accuracy 
from the plate solution. We excluded all 2MASS sources that do 
not have photometric quality flag from A to D in all bands. Also 
a few multiples and one clearly variable star were excluded and 
we were left with a total of 85 sources which could be compared 
in terms of fluxes and positions. The JHKs zeropoints were de- 
termined using 2MASS stars brighter than K s = 13.5 mag after 
excluding 2 cases of possible variability (i.e. 17 stars). The stan- 
dard deviation of the offset between NOTCam and 2MASS pho- 
tometry is 0.07, 0.05 and 0.05 mag in J, H, and Ks , respectively, 
for sources with Kg < 13.5 mag. We estimate the completeness 
limit of the observations to be roughly 18.5, 18.0, and 17.5 mag- 
nitudes in /, H, and K s , respectively. 

2.2.2. ISAAC H 2 observations 

Near-infrared observations, using the H2 filter centred at the v = 
1-0 line at 2.12yum, were conducted on 2000 November 10 using 
the ESO Antu (VLT Unit 1) telescope equipped with the short- 
wavelength arm (Hawaii Rockwell) of the ISAAC instrument. 
The ISAAC camera ( |Moorwood et al. |[T998] l contains a 1024 
x 1024 pixel near-infrared array and was used at a plate scale 
of 0.147 arcsec/pixel resulting in a field of view of 2.5 x 2.5 
arcmin 2 on the sky. At each of 5 different jitter positions (with 
the length of the jitter box equals to about half the field of view), 
one image with 60 s exposure time were taken. 

Standard procedures for near-infrared image reduction were 
applied (e.g. |Yun & Clemens | ([1994)) resulting in a final mo- 
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R.A. offset [arcsec] 

Fig. 1. Ha image towards IRAS 07141-0920 covering about 6 x 
6 arcmin 2 (5.4 x 5.4 pc 2 , for our estimated distance of 3.2 kpc). 
Notice the bright nebula and the patchy and wispy dark cloud 
material. 



saic image for each band. A sky image was computed for each 
dithered frame. Each frame was then sky subtracted and flat- 
fielded using master sky flats. The final clean images were cor- 
rected for bad pixels while constructing the final mosaic which 
covers about 5x5 arcmin 2 on the sky. Due to the presence 
of significant field distortion in these ISAAC images (two pix- 
els at the edges and 2.5 pixels in the corners), the individual 
frames were corrected before they were put together to form the 
mosaic. The correction for field distortion was performed us- 
ing IRAF/GEOTRAN together with the adequate correction files 
provided by ESO at their web page. The central, co-added re- 
gion, with enhanced signal-to-noise ratio covers about 2.5 x 2.5 
arcmin 2 on the sky. 



3. Extended nebula 

3. 1 . Morphology and colours 

Figure [T] shows the Ha image centered at the IRAS 07141-0920 
position and covering about 6x6 arcmin 2 . The high-resolution of 
ALFOSC allied to the very good seeing during the observation 
produced an image containing a wealth of detail. It reveals a 
large, extended nebula which appears to be illuminated by one 
or more objects possibly seen in this Ha image. 

Contrasting with the bright parts of the nebula, dusty, opaque 
patches are seen, specially in the northern and eastern sides. 
These dark clouds appear to be in the process of being shred- 
ded or blown apart by stellar winds, jets or outflows phenomena. 
They appear elongated, filamentary, wispy, clumpy or knotty. 
The overall image conveys an impression of streams, flows, tur- 
bulence, cavities, a region of great dynamical complexity. 

The main contributions to the nebula brightness are contin- 
uum ultraviolet scattered light (mostly contained in the U-band) 
and Ha emission. This is shown in the colour composite im- 
age shown in Fig. [2] This Figure shows a U-band (blue), R- 
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Fig. 2. U-band (blue), R-band (green), and Ha (red) colour com- 
posite image towards IRAS 07141-0920, covering about 6x6 
arcmin 2 . North is up and East to the left. Blue extended emission 
traces scattered light. Red extended emission traces Ha nebular 
emission from excited gas. The green colour (R-band contin- 
uum) appears essentially tracing point sources. 



band (green), and Ha (red) colour composite image towards 
IRAS 07141-0920. With this colour code, the blue extended 
emission traces scattered light while red extended emission 
traces Ha nebular emission from excited gas. The green colour 
(R-band continuum) appears essentially tracing point sources. 
Much of the nebular light seen is composed of Ha emission. 
The contribution of continuum scattered light is stronger near 
the centre of the nebula where the illuminator is most likely to 
reside. 
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B0.5V 




Fig. 3. Zoomed-in Ha image of the nebula covering 1 .9 x 1 .9 
arcmin 2 of the southern and central region of the nebula. B0.5V 
is the illuminator of the nebula. Notice the rim-bright features, 
some of them indicated by arrows. North is up and East to the 
left. 



by |Samal et al. | < f2007) >. In the remaining, we will refer to this 
object as S294-B0.5V. 

Additional arguments support this conclusion. Figure [5] re- 
veals nebular structures that display illuminated rims and shad- 
ows implying that the incident light is along the directions indi- 
cated by the two solid line arrows. In all cases, the only bright 
object located within the bright scattered light region is object 
S294-B0.5V. In section 5, we discuss optical and infrared pho- 
tometry results for this object, and the parameters of the associ- 
ated cluster. 



3.2. The illuminator star 

Given the presence of the HII region and the large, extended 
nebula, at least one O or B star is responsible for the ionization 
and excitation of the nebular gas. In order to identify the illu- 
minators), we note that it should be located close to the central 
region of the nebula which happens to be also where the blue 
scattered light is brightest. 

Because dark lanes of obscuring material can be seen in sil- 
houette against the bright nebula (mostly in the northern part of 
the nebula), we infer that the illuminator (or illuminators) must 
be partially embedded in the cloud. The illuminator should not 
be too deeply embedded, with its light escaping from "holes" 
(regions of lower optical depth, lower extinction) after scatter- 
ing off the walls of cavities in the near side of the cloud and 
exciting or ionizing the gas at the surface of the cloud. 

In this central region, the brightest source is the star marked 
with the letter "B0.5V" in Fig. [3] and we identify this object as 
being the illuminator of this nebula. This object coincides with 

T979) who first 



the star No. 4 of Moffat, Jackson, & Fitzgerald 



identified it as the hottest star in the region and thus the respon- 
sible for the HII region. This identification has been confirmed 



4. Stellar density distribution: a young stellar 
cluster 

Together with this massive B0.5V star, this region hosts the birth 
of a young stellar cluster associated to the massive star. 

Figure |4]presents the NOTCam JHKs colour composite im- 
age obtained towards IRAS 07141-0920. The optical nebula 
which dominated much of the UBVR and Ha images is absent. 
Instead, a large number of stars are seen, located in the region 
of the optical nebula but also to the east, where a dark patch of 
extinction is present in the optical images. 

A higher angular resolution image is provided by the VLT 
ISAAC Hi image seen in Fig. [5] The stars appear to form two 
distinct groups or sub-clusters: a western group (the western sub- 
cluster), that includes S294-B0.5V, less extincted and located to- 
wards the optical nebula; and an eastern group (the eastern sub- 
cluster), redder and located towards a dark optical region. 

The eastern part of Fig. [5] also shows the presence of elon- 
gated, jet-like or filamentary emission which suggest flows or 
material being blown away, likely by powerful stellar winds, or 
affected by ionization fronts. 

The presence of the two groups of stars is more clearly seen 
in the contour plot of the stellar surface density, shown in Fig. [6] 
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Fig. 4. J (blue), H (green), and Ks (red) colour composite image 
towards IRAS 07141-0920. The image covers about 3.2 x 3.4 
arcmin 4 (dashed-line box in Fig. 7). North is up and East to the 
left. The white cross indicate the position of the IRAS source 
((0,0) position of Fig. [T) 




i o ^7 

R.A. offset [arcmin] 



Fig. 6. Contour plot of the stellar surface density of the H2 im- 
age. Two peaks of the stellar surface density can be identified: 
the strongest, in the West, is the "western peak" ( where the il- 
luminator is located, see text); the second one, in the East, is 
the "eastern peak". The axes give the coordinates relative to the 
IRAS source. The lowest contour level is 9 sigma above the sky 
background stellar surface density. The highest contour level is 
about 14 times the value of the sky background stellar surface 
density. 




i. A. ijf'sel [n-cmi - i] 



Fig. 5. Central region of the H2 mosaic image towards 
IRAS 07141-0920. The axes give the coordinates relative to the 
IRAS source. The two zoomed-in views reveal elongated emis- 
sion in the eastern part, and the enhanced concentration of stars 
around S294-B0.5V. 



obtained from the H2 mosaic image Two peaks of the stellar sur- 
face density can be identified. The stronger peak is located in the 
central region of the optical nebula and we name it the western 
stellar density peak. It contains S294-B0.5V, the illuminator of 
the optical nebula, and a steep slope in the number of stars sur- 
rounding it. Thus, in a small region around the massive S294- 
B0.5V star, there is a large number of young lower-mass stars 
(see also the zoomed-in upper-right panel of Fig. The stel- 
lar surface density at this western peak is about 14 times higher 
than the value of the sky background stellar surface density. The 
second peak is located in the dark region to the east of the neb- 



ula and we name it the eastern stellar density peak. Together, the 
two peaks correspond to the presence of a young stellar cluster 
showing sub-structuring in the scale of about 2 pc, the size of 
the region. The position of the IRAS source, at offset (0, 0), is 
located between the two peaks with the IRAS uncertainty ellipse 
along the east-west direction. 

Interestingly, the map of 1.46 GHz continuum emission 
( Fich ||T993] > also displays two peaks located approximately at 
the same positions of the peaks of the stellar surface density. 
However, the stronger radio continuum peak is associated with 
the eastern peak, that is with the more embedded sub-cluster, in- 
visible in the optical images. The 1280 MHz radio continuum 
map of Samal et al. (2007 1 contains a peak coincident with the 
stronger 1 .46 map, and thus the eastern stellar density peak. 

5. Photometric results 

5.1. Standard calibration 

As mentioned above, observations of standards stars in the 



Landolt catalogue ( |Landolt 1992) were secured for standard op- 
tical calibration purposes. On the other hand, when additional re- 
liable standard values are available for a region of the sky, an ad- 
vantageous calibration procedure is at hand, by means of direct 
correlation between the instrumental values and those published 
standard va lues ( [Delgado & Alfaro ||2000l |Delgado, Alfaro, &| 



Yun 20071. The only published photometric study including 



stars of this region (Moffat, Jackson, & FitzgeraId~|([T979); MFJ 



in the following) contains UBV indices for 15 stars included 
in the field of our UBVRI observations. Of these 15 stars, we 
have good, non-saturated photometry for 14 stars. The correla- 
tion with the values from this study can be performed with some 
caution. Several stars included in this photoelectric study are lo- 
cated in the central region of the nebula, and their diaphragm 
photometry could be affected by the nebula brightness and some 



6 



Yun, Djupvik, Delgado, & Alfaro: A young double stellar cluster.. 



undetected crowding. In particular, this could be the case of the 
most luminous star, S294-B0.5V (MFJ No. 4), which is sur- 
rounded by several less luminous stars, more clearly visible in 
NIR frames. To avoid the influence of these effects, we used, for 
the correlation, 8 out of the 14 common stars, which are located 
outside the small central region marked in Fig. [7] This area en- 
closes the concentration of stars seen in the UBVRI frames, and 
more clearly in the NIR frames. The resulting correlation be- 
tween UBV indices is very well defined, and provides a set of 
standard, MFJ-derived, (U - B), (B - V) and V values. 

These values derived from photoelectric photometry are pre- 
sumably well tied to the original Johnson system. They are, 
however, systematically different from those obtained through 
calibration with Landolt standards. For the 8 stars used in our 
correlation, we obtain Vumdoit - Vmfj = -0.04 ± 0.02; (U - 
B)Landoi, ~(U- B) MF] = -0.07 ± 0.03; (B - V) Wo „ - (B - 
V)mfj - 0.01 + 0.03. These differences between UBV Johnson 
and Landolt systems were previously found and discussed by 
Delgado & Alfaro (2000). In our case here, the ZAMS line which 
we use for color excess and distance determination ( Schmidt - 
Kaler fT982| l is likely to be more adequate for Johnson values, 
and we adopt here the Johnson-derived set of standard indices. 
As for the calibration of (V - R) and (V - I) indices, it is ob- 
tained from the Landolt standard star measurements, following 

— < 1998). The 



IRAS_07141-0920 



the multilinear model described by Delgado et al. 



uncertainties of all the standard indices obtained, given by the 
root mean squared deviation of the residuals, amount to 0.02, 
0.03, 0.03, 0.06, and 0.03, for V, (U - B), (B - V), (V - R), and 
(V - I), respectively. 
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Fig. 7. Schematic representation of the field covered by the 
U BVRI observations. Only stars detected in the V-band are plot- 
ted. The smaller central square (solid-line) defines what we name 
core-stars (see text). The larger dashed-line box indicates the 
field of the JHK S NIR images. 



5.2. Colour excess, distance, and age 

The standard magnitudes of all the stars detected in our 
UBVRIHa images are listed in Table l|^] In this Table, columns 
(1), (2) and (3) identify the stars giving respectively, the run- 
ning ID number, and the equatorial coordinates (R.A. and Dec, 
epoch 2000); columns (4) through (14) give the photometric in- 
dices V, (U - B), (B - V), (V - R), (V - I), R, I, (R - Ha), 
J, H, and K$, and their corresponding uncertainties in paren- 
thesis. Column (15) gives the cluster membership assignments 
from optical photometry (see section 5.3). Column (16) indicates 
whether the source has near-IR excess. Finally, column (17) con- 
tains the identifications of the stars if it appears in the WEBDA 
data bas^] 

A schematic chart of the field covered by our optical obser- 
vations is plotted in Fig. [7] where all stars with V measurements 
are shown. In the following discussion, we will first concentrate 
on the stars located in the smaller square region at the center of 
the field, indicated in Fig. |7j For simplicity, we refer to these 
stars as the core-stars. 

In Fig. [8] we plot the Colour-Colour (CC) diagram (U-B) vs 
(B - V), and the Colour-Magnitude (CM) diagrams V vs (B - V) 
and V vs (V - /) of all the observed stars. Among these stars, 
the core-stars from Fig. 7 are represented by open circles and 
open squares. The latter appear to be part of a group of fore- 
ground stars at very low reddening. In the CC diagram (fig [8j, 
the ZAMS line (Schmid t-Kaler |19 82) shifted for a color excess 
of E(B-V)=0.1 (dashed line) fits well this group of stars. These 



Table 2. Parameters for nebula stars 



4 Table 1 is available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr or via http://cdsweb.u-strasbg.fr/cgi- 
bin/qcat?J/A+A/vol-number/page-number 

3 This research has made use of the WEBDA database, operated at 
the Institute for Astronomy of the University of Vienna. 
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foreground stars are homogeneously distributed in the observed 
field, and thus include some core-stars, namely, those we repre- 
sent by squares. 

On the other hand, the open circles appear to be a sequence 
of reddened core-stars. In fact, the solid line is the previous 
ZAMS line shifted to a location defined by E(B-V)a:1.3, with a 
reddening slope a=0.72. With these parameters the solid line fits 
well the shape of the most defining part of the sequence, namely 
that close to late B and early A type stars. Thus, we assume that 
the stars represented with open circles are very likely to be red- 
dened stars associated with the nebula. They correspond to the 
visible members of the young stellar cluster, better identified in 
the spatial distribution seen in the NIR images. We refer to them 
as nebula-associated stars, or just simply nebula stars. The two 
stars represented by crossed circles, which are indeed among the 
core-stars, are not considered in principle MS nebula stars, judg- 
ing from their locations in this diagram and in the CM diagrams. 
They belong to a conspicuous group of stars, which could be 
also associated with the nebula, and showing at the same time 
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Fig. 8. (U - B)-(B - V) CC diagram and V-(B - V), V-(V - I) 
CM diagrams. The so-called core stars (see text) are shown as 
squares and circles. The ZAMS line is plotted twice in the CC 
panel: with E(B - V) = 0.1 (dashed line, to account for the fore- 
ground field) and with E(B-V)=1.3 (solid line, to account for the 
sequence of reddened core-stars). The shifted line reproduces 
very well the sequence of reddened core-stars (circles), which 
we name nebula stars (see text). In both shifts the reddening 
slope E(U — B)/E(B - V) = 0.72 is applied. The two crossed 
circles are core stars rejected as MS nebula stars from their si- 
multaneous location in all three diagrams. The dotted line in the 
CC diagram encloses a region with stars possibly associated with 
the nebula and exhibiting (U-B) excess. Note also the distinctly 
deeper magnitude limits as we move towards redder bands, from 
U to /. 



a possible excess in the (U-B) index. This is a feature associ- 
ated with accretion in forming stars, as discussed by Gullbring 
et al. (1998, 2000). A more detailed discussion of these stars is 
postponed to the discussion of cluster membership. 

We are therefore left with 8 stars (listed in Table 2), which we 
consider as the most probable nebula stars. Table 2 also provides 
a quantitative estimate of the colour excess E(B - V) for the 
nebula stars, using the value of the reddening slope or=0.72. The 
calculated value for the nebula illuminator, star MFJ No. 4 (No. 
717 in Table 2), leads to a spectral type estimate of B0.5V for this 
star, in good agreement with MFJ and close to the spectroscopic 
type estimated by |Russeil, Adami, & Georgelin ( 2007| l. 

With the values of the colour excesses listed in Table 2, a 
determination of the visual absorption is possible, once the vi- 
sual absorption coefficient Ry = Ay/E(B - V) is obtained. Here 
we dispose of JHK S observations, which allow the estimate of 
Rv for individual stars with well defined E(B - V) color ex- 
cess. For a particular star, an E(V-K) value is calculated in the 
V-K vs B-V diagram, by shifting the ZAMS line (Koornneef 
H |1983| l in the B-V axis with the known E(B-V), and calculat- 
ing the distance to the line in the V-K axis. For our nebula stars 
we then obtain Rv values through the relation Rv(K)=1.12*E(V- 

These values provide an 



1999! 



K)/E(B-V)+0.02 ( |Fitzpatrick 
average Ry(K) = 4.5 + 0.6 from our NIR observations, and 
R V (K) = 4.3 ± 0.6 with K magnitudes from 2MASS. These high 
values of the absorption coefficient are however most probably 
affected by the NIR properties of the nebula stars. At least two of 
the nebula stars in Table 2 (520 and 660) show clear signs of NIR 
excess. In the presence of probable, or even suspected K$ -band 
excess, the use of the above relation is not adequate. On the other 
hand, a more reliable estimate of the absorption coefficient can 
be obtained from (V-I) photometry (Guetter & Turner ( 1997 ); 
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Fig. 9. (V - /) vs (B - V) Colour-Colour plot. Symbols for core- 
stars and nebula stars are as in Figure 8. The ZAMS line is 
shifted for a color excess of (E(B - V)=1.26, and a reddening 
slope E(V-I)/E(B-V)=l. 44. 



thors, R v (l) = 2.40 x 



E(V-I) 
E(B-V)- 



is valid for (V - I) indices in the 



GT97 in the following). The simple relation used by these au- 



Rron-Cousins system. The resulting Rv values are adopted here. 
In Fig. 9] we have plotted the CC diagram (V - /) vs (B - V). 
Symbols have the same meaning as in Fig. [8] The line plotted 
represents the ZAMS ( Schmi dt- Kaler |1982| l, shifted to account 
for the calculated average value of E(B - V), and a reddening 
slope J^gTm = 1.44. We can see that this value of the slope is 
adequate for all nebula stars. The resulting values of E(V - I) 
and Rv(I) are listed in Table 2. 

The values of the distance modulus obtained are also listed 
in Table 2. The average values E(B - V) = 1.26 ± 0.06, and 
DM=12.5 ±0.5 are obtained for the nebula stars listed in Table 2. 
Even with the large dispersion, this value is distinctly smaller 
than the distance given by MFJ, but agrees well with the re- 
vised value of 3.2 kpc obtained by Russeil, Adami, & Georgelin 
[1 ( |2007[ l. The larger distance obtained by MFJ is due to their 
adopted value for the absorption coefficient, around Ry-3.1, too 
small when compared to our calculated values. The presence of 
R v values both smaller and larger than those describing the av- 
erage galactic extinction have been reported in the literature (see 
[Turner] ([1994); Delgad o" et al. | (|T998) and references therein). 
GT97 find a similar value, and discuss their significance as rep- 
resentative of the extinction law in their line-of-sight. They ar- 
gue that the peculiar absorption coefficient indeed reflects the 
extinction law in their line of sight, rather than being related to 
the material physically associated to the stars. In our case, the 
dispersions of the averages for reddening slope and absorption 
coefficient are higher than theirs, partly because of larger errors 
in our photometry, partly because some of the peculiarities of the 
extinction in the region could be indeed related to specific prop- 
erties of the intracluster or circumstellar material in our region. 
This possiblity is suggested by the larger Ry(K) values found in 
our case than the ones by GT97, whereas the E(V - I) excesses 
and Ry(I) values are similar to theirs. 

We notice that, in general, the presence in a region of Ry 
values higher than the standard value (Ry=3.1) adopted for the 
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Fig. 10. CM diagrams V vs (B - V) and V vs (V-I), with indica- 
tion of cluster member candidates. Circles are MS and post-MS 
members. They include the nebula stars in Table 2. Large dots 
are PMS members assigned in the two CM diagrams V vs (V-R) 
and V vs (V - I). Squares are PMS members assigned in three 
CM diagrams, V vs (B - V), V vs (V - R) and V v s (V - I). The 
LogAge(yr)=6.6 (- 4 Myr) post-MS isochrone of Girardi et al. 
(2002 ) is plotted (solid line), shifted to account for the colour ex- 
cess and distance derive d for star S294B0.5V (MFJ No.4). The 
Palla & Staller | ([1999} isochrones for 1 and 10 Myr are also 
plotted (dashed lines), shifted for the average colour excess and 
distance modulus computed from the values of the nebula stars 
(see Table 2). Isochrones are labeled with their ages in Myr. 



average galactic extinction, might be caused by absorption due 
to material associated to the observed stars, and is usually men- 
tioned in connection with the observation of young star clusters 
( |Turner |1994)|de Winter et al. |19"97] >. On line with this reason- 
ing, the total reddening in the direction of the object would be 
composed of a foreground component which can be assumed to 
follow a normal reddening law, characterized by Ry = 3.1, and a 
component originating in the intracluster material, associated to 
the young stars, and responsible for the increase of the Ry value 



( |Pandey et al p003j|Sharma et al. |2007||Pandey et al p008] l. In 
the present case, the evidences from our CC diagrams (Figs. 8 
and 9) indicate that only a small part (0. 1 to 0.2 dex in E(B - V)) 
of the reddening affecting the nebula stars is to be considered as 
due to foreground reddening. Assuming that the absorption due 
to this part is described by Ry=3.1, would have little influence 
on the final distance estimate. 



5.3. Cluster membership: optical photometry assignments 



Using the p rocedures described in detail by Delgado, Alfaro, &| 
Yun ( 2007 1, we have analysed whether each star is likely to be- 
long to the young cluster. This membership analysis is applied to 
the optical UBVRI observations. It consists of calculating colour 
excesses and distance moduli for all stars with respect to the 
ZAMS and PMS isochrones. A star is then considered a clus- 
ter member candidate (MS or PMS), if its colour excess, and 
distance modulus, with respect to any of the comparison lines, 
coincide (within the errors) with the corresponding average val- 
ues calculated for a previously selected number of bona fide un- 



evolved MS members. In the present calculation, the MS mem- 
bers are the 8 nebula stars mentioned above, which provide the 
values listed in Table 2, and the corresponding average value of 
colour excess and distance modulus of the cluster. 

The reference lines used in the comparison are the ZAMS 
line (Schmidt- Kaler 1982|l, and th e PMS isochrones from two 
different models ([Pal la & Staller | fi99"9"l >, and |Siess, Dufour, &| 



Forestini (2000), in the following P99 and S00, respectively) 
When the distance modulus for a star is calculated from com- 
parison with PMS isochrones, we obtain a PMS membership 
guess for it. The comparison can be performed in several CM 
diagrams. Some differences occur when using P99 or S00 mod- 
els: the number of stars that fullfil the membership criteria are 
different: 58 stars turn out to be candidate PMS members with 
respect to P99 (assignment in 2 or 3 CM diagrams), whereas only 
38 fullfil the requirements when compared to S00. The bulk of 
this difference arises in the number of stars assigned as mem- 
bers around the spectral types A and F. Assigned PMS members 
of spectral types around G and later coincide for both isochrone 
models. The difference is due to the fact that the S00 isochrones 
are more luminous (around 0.5 brighter in My around A and F 
spectral types) than the P99 isochrones. As a result, also the av- 
erage age of the PMS candidate members is younger when mea- 
sured using the P99 isochrones. On the other hand, practically 
all stars assigned with respect to S00, are also assigned with re- 
spect to P99 models. The assignments listed in Table 1 were con- 
sidered with respect to either of the two models. The entries in 
column (16) are: 0, for no assignment; 1, for MS members; 2, 
for PMS members assigned in 2 CM diagrams; and 3, for PMS 
members assigned in 3 CM diagrams. 

The CM diagrams in Fig. 10 show the membership assign- 
In these plots, we show the 
[2002> , 



ments obtained in the analysis. 
post-MS isochrone of logAge(yr)= 



=6.6 ( Girardi et al. 
shifted to account for the parameters of color excesses derived 
for S294B0.5V. PMS isochrones of 1 and 10 Myr from P99, 
shifted to account for the average parameters derived above, are 
also plotted. Stars indicated with larger dots in Fig. 10 are those 



with a PMS membership assignment in two CM diagrams (V 
vs (V - R) and (V - I)). This means ( |Delgado, Alfaro, & Yun~| 
2007 ) that both colour excess and distance modulus calculated 
with respect to some PMS isochrone between 1 and 10 Myr, 
coincide (within errors) with the average values obtained from 
the nebula stars. Similarly, stars assigned as members in three 
CM diagrams are marked as squares. We detect PMS candidates 
down to masses around 0.7 M . The average age value for these 
PMS candidates can be taken as an age estimate for the possible 
PMS sequence. The values obtained for stars with membership 
assignments in at least three CM diagrams, are 7.4+1.3 Myr for 
P99 isochrones, and 8.3+0.7 Myr for S00 isochrones. These val- 
ues are higher than the age estimate obtained for the MS mem- 
bers. The significance of this estimate has yet to be assessed. If 
confirmed, it would indicate a sequential star formation in the 
region, with the most massive stars being formed last. This kind 
of trend is also suggested by the results on other regions of active 

2002]l ~" 



formation ([Mamajek, Meyer, & Liebert 
& Dias | ( [2003| l, |Delgado, Alfaro, & Yun \p.QQ7\ ) 



Sartori, Lepine. 



5.4. Cluster members with near-IR excess 

The J-H/H-Ks diagram for 437 NOTCam sources with cr Ks < 
0.2 mag is shown in Fig. [TT| The loci o f giant, supergiant, and 
main-sequence stars ( Koornneef 1 1 983] > are indicated with bold 
curves. We have used the Ax oc A parametrization of the NIR 
extinction law (Whittet 1988) to calculate a reddening slope of 



Yun, Djupvik, Delgado, & Alfaro: A young double stellar cluster.. 



9 




Fig. 11. J - H/H - Ks diagram for the 437 sources with cr Ks < 
0.2 mag. The reddening slope of 1 .6 results from applying the 
A,i oc /T 1 7 parametrization of the NIR extinction law and the 
reddening vector of an AO star is drawn (dashed line). The bona- 
fide IR-excess sources (big dots) are located more than 2 cr to 
the right of and below the reddening line and the probable IR- 
excess sources (triangles) have a separation from the reddening 
line between 1 and 2 cr of the individual errors in the J - H and 
H-K$ indices. Large circles mark nebula stars listed in Table 2. 



1.6 in the J - H/H - K s diagram for the NOTCam JHK S filter 
pass bands at AA 1.247, 1.632, and 2.140 //m, respectively. The 
reddening vector for an AO star is plotted as a dashed line. Most 
stars are located around the reddened loci of main-sequence and 
giant stars. 

Sources located right- and downwards of the reddening vec- 
tor of an AO star have excess emission in the infrared, a clear 
sign of hot circumstellar dust. We define as bona-fide IR-excess 
sources those separated from the reddening vector down and 
right by more than 2cr of the individual errors in the {J - H) and 
(H-Ks) indices (48 sources) and as probable IR-excess sources 
those with a separation between l<x and 2<x (17 sources). These 
are marked as big dots and triangles, respectively, in Fig. 1 1 The 
8 encircled sources in the figure are the so-called nebula stars 
discussed in Sect. 5.2. Two of them have clear or probable near- 
IR excesses (IDs 520, 660) and one (ID 733) is uncertain (see 
Table 2). In Table 1 the column IREX has the value for no IR- 
excess, 1 for probable IR-excess and 2 for bona-fide IR-excess 
sources. The sources with excess emission in the infrared ow- 
ing to hot circumstellar dust are candidate Young Stellar Objects 
(YSOs), either pre-main sequence (PMS) stars with circumstel- 
lar discs or protostars embedded in a more spherical dust distri- 
bution. We cannot separate clearly between PMS stars and proto- 
stars only from near-IR photometry, but because of the different 
time scales of these populations, the majority of the IR-excess 
YSOs are probably PMS stars with discs. 

There are six very red objects near the centre of the embed- 
ded cluster which are not detected in the J band. The lower limit 
(J - H) values suggest that three of these have IR excesses and 
are probably very embedded YSOs in their earliest evolution- 
ary phase. These have been given the IREX indicator value 1 in 
Table 1. 

Only 144 (272) of the 444 near-IR sources have V {I) 
band counterparts, and these are typically the less extinguished 
sources with Ay less than 10 mag approximately. Of the 68 IR- 
excess selected YSOs only 10 (30) are detected in V (I), and 
of these, 2 have been designated PMS stars from the optical 
photometry fulfilling the critera of two CM diagrams (ID 412, 



473, 822). Also, 3 of the 12 optically selected MS stars have 
IR-excesses (ID 343, 520, 660). We note that our near-IR obser- 
vations cover only the central 3.7' x 3.4' part of the 6.5' x 6.5' 
optical field (see Fig. 7). 

The analysis of the optical photometry finds a population of 
58 PMS sources, of which only 2 overlap with the 68 IR-excess 
YSOs. This shows that the two photometric methods are comple- 
mentary and sample different parts of the YSO population of this 
type of regions. The optically selected PMS stars are in general 
expected to be closer to the main-sequence phase and have less 
or no measureable circumstellar dust. It is also known that only 
about 60% of the Classical T Tauri stars in Taurus-Auriga have 
measurable NIR excesses (see Fig 7 in Strom, Strom, & Merrill 
|1993| l. In addition, a large fraction of the IR-excess population 
is not optically visible, making a direct comparison difficult. 

The fact that 25% of the optically selected main-sequence 
stars show IR-excesses is probably related to the fact that these 
are all relatively massive, and in theory they arrive on the main- 
sequence already while still accreting material, skipping the 
PMS phase. It is also possible that these 3 out of 12 main- 
sequence stars are unresolved binaries. 

The IR-excess population may not be well sampled by near- 
IR photometry, however. The mid-IR ISC^Jsurveys of young em- 
bedded clusters demonstrated well that in all regions only about 
50% of the IR-excess population is sampled by near-IR photom- 
etry (Kaas & Bontemps 2001). Mid-IR photometry can better 
separate IR-excess from reddening and detect hot dust which is 
obscured from view. 



5.5. The luminous protostar candidate 

Searching the IRSA archive we found one mid-IR source in- 
side our near-IR NOTCam field in the MSX6C catalogue from 
the Midcourse Space Experiment. This is G224. 1880+01.2407 
which is located within 1 .6" (with position angle 255 degrees) 
from a IR-excess source in our sample (ID = 1050, Ks = 15.1 
mag). This is less than 3cr of the MSX6C positional uncertainty, 
and therefore we consider the source ID=1050 a possible near- 
IR counterpart. According to its Spectral Energy Distribution 
(SED) in the mid-IR (see Fig[T2|), it is a protostar candidate. The 
SED inde^ from 2 to 15 fim is as high as 0-2-15 = 2.7, which ac- 
cording to the YSO classification scheme (Lad a |1987| Greene 



|et al. 1 1994] ), though strictly valid for low mass stars only, means 
this is a protostar candidate of the type Class I or possibly in 
the transition phase from Class to Class I. The brightness of 
the object in the mid-IR (0.8755 Jy at 8.28 //m) indicates it is a 
luminous/massive protostar. 

Combining the NOTCam JHKs fluxes with the fluxes in the 
MSX bands ACDE at AA 8.28, 12.13, 14.65 and 21.34 /mi, re- 
spectively, and the IRAS 12, 25, 60 and 100 //m fluxes, the SED 
is sampled over a sufficient wavelength interval to compare it 
to YSO models (see Fig. [12). The online YSO SED model fit- 
ting tool provided by Robitaille et al. (2007 ) was used to find 
the best YSO model for the observed SED. The input ranges 
in interstellar extinction and distance was limited to 10 mag 
< Ay < 40 mag and 2.1 kpc < D < 4.1 kpc, respectively. The 
models which best fit the observed SED all give: a stellar mass 
from 8 to 12 M Q , a total luminosity of 1 to 4 10 3 L Q , a stel- 
lar age of 40000 yrs or less, and an envelope accretion rate of 



6 Infrared Space Observatory 

7 The SED index is defined as ff2-i5 



log(v 15 F v )-log(v2F V2 ) 



where F v 



l0gVl5-l0gV2 

are the flux densities (in Jy) at the frequencies v (in Hz) corresponding 
to the wavelengths 2 and 15 jim. 
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Fig. 12. The Spectral Energy Distribution (SED) of the luminous 
protostar candidate sampled by NOTCam JHK$ band fluxes, 
MSX ACDE band fluxes, and IRAS 12, 25, 60, and 100 /mi 
fluxes. The best fitting model found using the online SED fit- 
ter ( Robitaille et al. |[2007[ ) is plotted as a continuous line that 
splits in two at longer wavelengths for small (lower) and large 
(upper) aperture sizes. The Y-axis is scaled to a distance of 1 
kpc. 



~ 5 x 10~ 3 M Q /yr. Adjustments for aperture mismatches of the 
different wavelengths do not significantly alter the derived val- 
ues. 

For all best fitting models, the measured flux in the MSX D 
band (13.5-15.9 /mi) is too low. We speculate that this is caused 
by very strong absorptions, such as the known feature at 15.2 
fim, found by Alexander at al. ( 2003| > to be present for the 
youngest and most embedded YSOs they surveyed with the 
ISOCAM CVF. The feature is believed to be attributed to the 
bending mode of CO2 ice. 

An MSX A band image (6.8 - 10.8 /mi) covering the area 
of our optical and near-IR studies is shown in Fig. 13 The error 
ellipse of the IRAS point source is shown, and for reference, the 
position of the massive main-sequence star S294-B0.5V (the op- 
tical nebula illuminator) is marked with a small circle. We note 
that in addition to the bright MSX point source, the 8/mi image 
shows large scale diffuse emission extending as curved arms ~ 2' 
or more from the point source. The diffuse emission is strong at 
8 //m, but barely distinguished from the noise in the other MSX 
bands. Its main origin is likely the emission features at 7.7 and 
8.6 /mi, two of the emission bands of polycyclic aromatic hy- 



drocarbons (PAHs). As shown in Fig. 15 the curved filaments in 



emission at 8 /mi trace well the borders of the Ha emission. This 



resembles very much the findings of Deharveng et al. (2003) for 
the geometrically simpler sources Sh 217 and Sh 219, two spher- 
ically symmetric HII regions surrounded by annular photodisso- 
ciation regions traced by both 21cm HI and 8//m PAH emission. 
In our case the morphology is more complex. 

5.6. PMS membership and Ha emission 

The combined evidence from all photometric diagrams shows 
indications of several stars associated to the cluster, which are 
not detectable as candidate members at optical wavelengths. We 
consider here spectral features which could be present in several 
stars in our sample, and have been identified in the literature as 
spectral signatures characteristic for PMS stars. The basic fea- 
ture in this context is the presence of Ha in emission, typical in 



cr 
7-17- 


I0.p 50 


, , ,40 


» i i i 30 


, , ,20 


• 16:1 


o-q— 


24:0 




• m 










0.0 26:0 














0.0 28:0 














0.0 30:0 






























Fig. 13. MSX A band (8.28 /mi) image with the error ellipse of 
IRAS07 14 1-0920 taken from the IRAS point source catalogue 
as well as the location of the optical star S294-B0.5V (circle). 
The FOV is roughly as for the optical images, and the J2000 RA 
and DEC grid is shown. The intensity is in logarithmic scale. 



the spectra of low mass forming stars of the TTauri type, both 
in classical (CTT) and weak-lined (WTT) stars, as well as in 
higher mass forming stars of the HAeBe type. The value of the 
R - Ha index has been used in previous studies of young clus- 
ters to describe this feature, although its value may be reflecting 
spectral properties other than those due to the PMS nature of 
t he star s ([Sung, Bessell & Lee | ( |1997| ), |Sung, Chun, & Bessel| 
[] ( |2000[ ), |Sharma et al. | ( |2007| l). For our stars this index, cal- 
culated as the difference between instrumental magnitudes, is 
listed in column 11 of Table 1, and plotted versus V magnitude 
in Fig. 14 Different symbols are used to distinguish stars which 



show different features or spectral signatures of interest: possible 
Ha emission (R - Ha > -3.4 in the plot), PMS candidates from 
our UBVRI photometry, NIR excess, as described in Section 5.4, 
and probable excess in U - B (see Sect 5.2 and Fig. 8). 

In their spectroscopic study of the PMS population in the 
4 Myr old cluster Tr 37, |Sicilia-Aguilar et al. ] ( p005] l find 40% 
of stars with signs of accretion, a percentage that decreases to 
almost zero in the 10 Myr old cluster NGC 7160. Furthermore, 
the percentage of these accreting stars with signs of NIR excess, 
50% in Tr 37, drastically decreases with increasing age. The cor- 
responding numbers from our photometric analysis would be: 16 
stars with (R - Ha) > -3.4 out of 58 stars selected as optical 
PMS candidates (27%), and 4 stars among these 16 (25%) with 
NIR excess. Considering Ha emission as a possible tracer of 
accretion processes, and given the age of our optically selected 
PMS sequence (around 8 Myr) the corresponding percentages in 
our sample are consistent with those results. 



Figure 14 indicate that the presence of possible (U - B) ex- 
cess coincides frequently with the indication of Ha emission, 
but it seems to exclude PMS candidates selected in 2 or 3 optical 
CM diagrams. In particular, those of earlier spectral type (around 
AF: stars marked with both an open and a closed circle), most 
do not show any sign of circumstellar activity or ongoing ac- 
cretion, whereas several fainter PMS candidates, below spectral 
type G and later (marked with only open circles), do show signs 
of Ha emission, but only two of them possible (U - B) excess. 
This mutual exclusion of PMS nature (inferred from optical CM 
diagrams) and (U - B) excess might have two possible causes. 
Stars with (U-B) excesses might have other spectral anomalies, 
such as (B - V) indices with some excess, and possible Ha emis- 
sion. These features would prevent them from being selected as 
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Fig. 14. Plot of the instrumental index (R - Ha) versus V mag- 
nitude. Small dots represent all stars in our sample with valid R 
and Ha measurements. Different symbols are used as follows: 
squares, for stars with NIR excess; crosses, for stars with possi- 
ble emission in Ha ((R - Ha) > -3.4); open circles, for PMS 
candidate members, with respect to P99 or S00 models, selected 
in 2 CM diagrams; closed circles, PMS candidate members se- 
lected in 3 CM diagrams; triangles, for stars with possible (U-B) 
excess. 



optical PMS candidates by our CM diagrams method. On the 
other hand, and more interestingly, the disappearance of ( U - B) 
excesses as we move to earlier type stars could be due to mag- 
netospheric accretion stopping operating in HAeBe stars due to 
weaker or absent magnetic fields (Wad e et al. |2007[ ) when com- 
pared to magnetic fields in lower-mass PMS stars. 

In general, the progressive disappearance in our PMS sample 
of spectral features, as Ha emission, and even (U-B) excess, 
as we move to earlier spectral types, is a finding on line with ar- 
guments in favour of a faster dissipation of accreting discs and 
envelopes in stars of higher masses (Herna ndez et al. |2005] >. In 
fact, a decrease in number and intensity of emission features in 
stars of spectral type F and earlier has been found in compre- 
hensive analyses of the spectral types of PMS stars ( Hernandez 
|et al. ||2004| ). Furthermore, the results of |Manoj et al. | ( |2006] >7 
relative to HAeBe stars, show a rapid decrease of Ha emission 
with age, but more interestingly, no F-type star in their sample 
exhibits an equivalent width of Ha in emission above 10, inde- 
pendently of age. In this same context, Hernandez et al. ( 2006 ) 
have found that the PMS population in Orion, between ages of 5 
and 10 Myr, show signs of residual debris discs, without excess 
in the IRAC bands and without emission lines in their optical 
spectra, but with a varying degree of 24 yum excess. 



6. Discussion 

6.1. Spatial distribution 

A total of 133 cluster member candidates were found in 
the previous section using CM and CC diagrams based on 
UBVRIJHKs photometry. From the optical CM diagrams us- 
ing the method of Delgado, Alfaro, & Yun ( |2007| l we find 12 
main-sequence stars and 58 PMS stars. Using IR-excess as a 
youth and member criterium we find 68 PMS members from 
the J - H/H - Ks diagram. Only 5 of these coincide with the 




Fig. 15. The spatial distribution of the cluster members overlaid 
on the Ha image (grayscale). The white contour lines repre- 
sent the MSX 8 //m isophote emission. The positions of main- 
sequence stars (blue diamonds), PMS stars selected from optical 
CM diagrams (yellow circles), IR-excess stars selected from the 
J - H/H - Ks diagram (red squares), and Ha emission sources 
selected from the (R-Ha) index (green crosses) are shown. The 
massive MSX protostar candidate is located at the strongest peak 
of the 8 /mi emission. The spatial distribution of the IR-excess 
sources is strongly sub-clustered. 



optically selected members. In addition, we find 43 sources with 
strong Ha emission using the R-Ha index. Of these, 21 coincide 
with members selected from optical CM diagrams or with NIR 
excess. The total number of cluster member candidates found in 
this study is thus 155. 

The spatial distribution of the cluster members can give valu- 
able information about cluster evolution and star formation sce- 
nario. In Fig. [15] we show the location of the cluster member 
candidates with different symbols according to selection method 
(see figure caption for details) overlaid on the Ha image. While 
it is clear from Fig. [6] that there are mainly two stellar density 
peaks at 2 /mi, referred to as the eastern and western clusters, 
the spatial distribution of the cluster member candidates, as se- 
lected by our previously described methods, gives a more nu- 
anced impression. The density of cluster members is still highest 
at the location of the eastern and western cluster, but about half 
of the members are extended both towards SE and NW, tracing 
the dense middle lane of the region seen in absorption on the 
Ha image. On a large scale, star formation seems to be going on 
along a filament or flattened density structure along a SE-NW 
oriented ridge. 

As many as 8 of the 12 optically selected MS cluster mem- 
bers are located in the area of the strongest (western) stellar den- 
sity peak seen at 2.12 /mi (see Fig. |6]l around the source S294- 
B0.5V. Three are located in an area about 1.5' to the NW, and 
the last one about 1 .5' to the NE. We note that there is only one 
MS member found towards the SE. 

The optically selected PMS cluster members are mainly dis- 
tributed in roughly three loosely clustered areas: one around the 
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exciting source S294-B0.5V, one about 1.5' to the NW, and one 
about 1.5' to the SE. 

The IR-excess selected YSOs have a spatial distribution 
mainly across the NW-SE diagonal of the 3.7'x3.4' JHK$ image. 
In addition to the more scattered component of IR-excess source 
distribution, about half of them are strongly clustered. Around 
the source S294-B0.5V and coinciding with the strongest stel- 
lar density peak of Fig. [6] as many as 8 IR-excess sources are 
inside ~ 16" (see the zoom-up in Fig.|5]for high spatial resolu- 
tion). Towards the eastern cluster, there is clustering of IR-excess 
sources in a dense part of the cloud practically void of optical 
sources. The projected separations between all the near-IR ex- 
cess sources have a distribution with one peak at about 15 arcsec 
and another at 90 arcsec. The latter reflects the extent of the most 
active region of star formation, about 1 .4 pc for an adopted dis- 
tance of 3.2 kpc. The peak at 15 arcsec or 0.23 pc reflects the 
typical size of the smallest stellar groups. This is in good agree- 
ment with the smallest clustering scale of Class II sources (0.25 



pc) fo und in the nearby Serpens star forming region (Kaas et al. 
2004). In the easternmost of these groups resides the luminous 
MSX source G224. 1880+0 1.2407, a massive protostellar candi- 
date, the youngest of the members. 

Whether the groups of the embedded eastern cluster are 
younger than the members in the western cluster around S294- 
B0.5V is not clear. The clustering of IR-excess sources is strong 
around S294-B0.5 V as well. The age difference between the two 
most massive stars is clear, however. One is already on the main 
sequence, while the other one is still in its protostellar phase. It is 
also quite likely that the median age of the IR excess sources is 
less than the median age of the optically selected PMS sources, 
and for comparison the Serpens Class II sources with a similar 



clustering scale have an age estimate of about 2 Myr (Kaas et al. 
2004). Kinematic data may be needed in order to better under- 
stand the nature of the sub-clustering processes observed in this 
star forming region. 

6.2. A possible star formation scenario for Sh2-294 

The MSX A-band image shows diffuse emission along curved 
filaments that trace well the borders of the Ha emission (see 
Fig. 15 i. If the 8 /urn emission is mainly due to PAHs, as sug- 



gested in Sect. 5.5., this spatial distribution can be interpreted 
as a photodissociation region (PDR) forming a shell around the 
optical H II region, in a similar fashion as for Sh2-217 and Sh2- 
219 ( |Deharveng et al. | 2003 ), although the morphology is less 
symmetrical for Sh2-294. The bright MSX point source is ob- 
served in the direction of the eastern edge of the PDR region, 
only Aa = 24.3" and A5 = 4" offset from the eastern and 



strongest of the two peaks in the 1.46 GHz VLA map (Fich 
1993 ). The IRAS colours indicate an ultra-compact H II region, 
and the quoted 7' diameter of the VLA source ( Fich~| [T993" i is 
understood to include the free-free emission from the larger and 
optically visible H II region centred on the western radio peak 
close to S294-B0.5V. 

The appearance in the near-IR images is that of a double 
cluster with the eastern cluster containing the MSX protostar, 
probably responsible for the ultra-compact H II region, and the 
western cluster centred on the massive star S294-B0.5 V, respon- 
sible for the optical H II region. The fact that the mid-IR emis- 
sion traces so well the borders of the optical HII region strongly 
argues in favour of a physical connection and against a chance 
alignment of two clusters at different distances. Their projected 
separation is small, but the eastern, embedded cluster must be to- 
wards the front side of the optical H II region. Its shape formed 



by the sub-clustered IR-excess sources traces the east-west fil- 
ament (8/mi emission and Ha absorption). The spatial distribu- 
tion of the bulk of the cluster members, the lower mass stars, is 
not distinctly bimodal, however, and the optically selected PMS 
stars (with typical ages of 7-8 Myr) are in many areas spatially 
mixed with the IR-excess sources (most of which are not seen 
in the optical, are sub-clustered, and probably much younger). 
This can be partly due to projection and spatially variable cloud 
extinction. 

The on-going massive star formation, traced by the MSX 
protostar G224. 1880+01. 2407, seems to take place in the PDR. 
It may have been triggered by the optical H II region excited by 
S294-B0.5V, estimated to have spent about 4 Myr on the main- 
sequence. If this is the case, the age difference between the two 
generations is probably quite small. For the bulk of the cluster 
members, the low mass stars, it is hard to argue for any distinct 
age difference between the eastern and western part, although 
there is probably a difference between the median age of the IR 
excess sources and that of the optically selected PMS stars. In 
order to map the molecular cloud and get more insight into the 
star formation scenario for Sh2-294, follow-up observations in 
the mm and radio wavelengths are needed. 



7. Summary 

We present an optical and near-infrared imaging and photomet- 
ric study of the HII region Sh2-294 (S294) associated with the 
IRAS source PSC 07141-0920. 

High-resolution images in the UBVRI and Ha bands reveal 
an emission nebula with very rich morphological details. The 
nebula is composed of ultraviolet sctattered light and of Ha 
emission with local morphological features conveying an im- 
pression of a complex web of streams, flows, shocks, turbulence, 
cavities, a region of great dynamical complexity. Contrasting 
with the bright parts of the nebula, opaque patches are seen. 
These dark cloud regions appear to be in the process of being 
shredded or blown apart by stellar winds, jets or outflows phe- 
nomena. They appear elongated, filamentary, wispy, clumpy or 
knotty. 

Our optical photometry confirms that the illuminator of the 
nebula is likely to be star No. 4 of Moffat, Jackson, & Fitzgerald 
| (1979 1, with a spectral type of B0.5V, but located at a dis- 
tance of about 3.2 kpc, as recently revised by Russeil, Adami, 
&Georgelin1p007) . 

Our high-resolution near-IR images reveal an embedded 
cluster, extending for about 2 pc and exhibiting sub-clustering: 
a denser, more condensed sub-cluster (western sub-cluster) sur- 
rounding the high-mass B0.5V illuminator star; and a more 
embedded sub-cluster (eastern sub-cluster) located towards the 
eastern dark part of the nebula and including the luminous MSX 
source G224. 1880+0 1.2407, a massive protostellar candidate 
that could be the origin of jets and extended features seen at 2. 12 
/j.m. 

This double cluster appears to be clearing the remaining 
molecular material of the parent cloud, creating patches of lower 
extinction and allowing some of the least reddened members to 
be detected in the optical images. 

155 stars are considered as cluster members, 12 of them, 
MS members. The remaining are PMS candidates, inferred from 
three different types of indicators: membership assignment from 
comparison with PMS isochrones, Ha emission deduced from 
the (R - Ha) index, and presence of NIR excess. About half 
of the IR-excess sources are spatially distributed in sub-clusters 
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with typical sizes of about 15 arcsec or 0.23 pc. These are prob- 
ably younger than the optically selected PMS stars. The age es- 
timate for the massive MS member is about 4 Myr and the aver- 
age age of the optically selected PMS members around 7-8 Myr. 
Our results support time scale estimates for the duration of ac- 
cretion disks around 10 Myr for late G-type stars and later, and 
a remarkable decrease of this duration towards earlier spectral 
types. 

We provide a catalog containing the photometry of all 1344 
stars detected at least in one of the UBVRIJHKs bands. Cluster 
membership is also included. 
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